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1.  Introduction. 

If  the  characteristic  function  of  an  n-dimensional  random  vector  x 
has  the  form  exp(it,u)<j)(t,It) ,  where  y:  nx  1,  E:  n*n,  and  E  0, 
we  say  that  x  is  distributed  according  to  an  elliptically  contoured 
distribution  with  parameters  y,  E,  and  4>,  and  we  write  x  'v  EC  (y,E,<(>) . 

^  <v  -v  II  -v  <*• 

The  class  of  elliptically  contoured  distributions  has  been  studied 
by  several  authors:  Schoenberg  (1938),  Kelker  (1970),  Devlin,  Gnanadesikan 
and  Keltenring  (1976),  Kariya  and  Eaton  (1977),  Muirhead  (1980),  Cambanis, 
Huang  and  Simons  (1981),  and  Anderson  and  Fang  (1982). 

Statisticians  have  been  trying  to  extend  the  sample  theory  in  multi¬ 
variate  analysis  to  the  case  of  samples  being  dependent  or  the  case  of 
samples  being  from  nonnormal  populations.  In  this  paper  we  consider 
sampling  theory  in  which  the  distribution  of  the  population  belongs  to  the 
class  of  elliptically  contoured  distributions  and  the  samples  are  dependent 
According  to  this  requirement  multivariate  elliptically  contoured  distri¬ 
butions  are  defined  and  some  basic  properties  are  discussed  in  Section  1 
and  Section  2.  The  distributions  of  some  important  statistics  in  the 


sampling  theory  (such  as  the  correlation  coefficients,  the  multiple 
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correlation  coefficients,  Hotelling-T  ,  the  sample  covariance  matrix, 
the  generalized  variance,  the  quadratic  forms,  etc.)  are  obtained  in 
Sections  3  to  Section  6.  As  applications  of  the  theory  on  multivariate 
elliptically  contoured  distributions  we  consider  the  model  of  multiple 
regression  with  the  error  matrix  being  distributed  according  to  a 
multivariate  elliptically  contoured  distribution. 

Kariya  and  Eaton  (1977),  and  Muirhead  (1980)  discussed  the  effects 
of  elliptical  distributions  on  some  standard  procedures  involving  correla¬ 
tion  coefficients,  but  the  model  that  we  consider  in  this  paper  is  different 
from  theirs. 

Throughout  the  paper,  N^fti,!)  denotes  the  n-dimensional  normal 

2 

distribution  with  mean  y  and  covariance  matrix  £ ;  denotes  the  chi- 

squared  variable  with  k  degrees  of  freedom;  Bto^.o^)  denotes  the  Beta 

distribution  with  parameters  a,  and  a„;  D  (a, , . . . ,a  , ;a  )  denotes  the 

i  L  m  X  m-1  m 

Dirichlet  distribution  with  parameters  a, , . . . ,a  ;  F(k,£)  denotes  F-distri- 

J.  m 

bution  with  k  and  l  degrees  of  freedom,  tR  denotes  the  t-distribution 

with  n  degrees  of  freedom;  W^d.n)  denotes  the  Wishart  distribution  with 

covariance  matrix  E:  p*p  and  n  degrees  of  freedom;  U  denotes 

~  P,m,n 

Wilks'  statistic  which  is  the  ratio  J G {  / 1 G-l-H j ,  where  G  ^  Wp(E,n), 

H  w  (Z,m),  and  H  and  G  are  independent;  I  denotes  the  nxn 
P  ~  ~  ~n 

identity  matrix;  denotes  the  n*l  vector  with  elements  1;  rk(A) 

denotes  the  rank  of  the  matrix  A  and  A  denotes  a  generalized  inverse 
of  A. 


2.  Definitions  and  Basic  Properties. 

Let  X,  M  and  T  be  n*p  matrices.  We  express  them  in  terms 
of  elements,  columns,  and  rows  as 
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X  -  (x  )  -  (x. ,x„ , . .  • ,x  ) 
~  lj  ~p 


5(1) 


l*Cn)J 


,  x  -  vec  X'  , 


m  -  <w  )  - 


tyu1 


,  p  >  vec  M'  , 


I  “  (t1j)  "  (3l»52****»5p) 


3(1) 


~(n) 


,  t  -  vec  T' 


Here  x  -  vec  X'  ■  ^x(i)  ,x(2)  *  *  *  *  ,x(n)  ^  '  an<*  wit^  t*ie  same  meaning  for 
p  and  t. 


Definition  2.1.  If  the  characteristic  function  of  a  random  matrix 
X  has  the  form 

(2.1)  exp(l  ^^tJj)p(j))«|>(tJ1)Z1t(1),...,tJn)Znt(n))  , 

with  £-,...,£  >0,  we  say  that  X  is  distributed  according  to  a 

multivariate  (rows)  clliptlcally  contoured  distribution  and  write 
X  *v»  MK:  (M;!. 

—  u*p  —  -i 

Obviously,  whan  n  ■  1  the  multivariate  elliptically  contoured 
distribution  reduces  to  the  common  elliptically  contoured  distribution. 

Let  u(q)  denote  a  random  vector  which  is  uniformly  distributed  on 
tba  unit  skiers  in  ***  and  8^(1  t|p)  denote  its  characteristic  function. 


3 


Let  9  be  the  class  of  all  functions  $:[0,°°)  *  [0,®)  x  •••  x  [0,®) 

n«  t  •  •  • 

1  n 

2  2 

■  [0,°°)n  +  3R  such  that  <K|ltjJ|  ,...,||t  ||  )  is  a  characteristic  function, 

where  t^,...,tn  are  x  1, ... ,mn x  1  vectors,  respectively. 

By  a  method  similar  to  Schoenberg  (1938) ,  it  can  be  shown  that 

$  €  9  _  if  and  only  if 

®-i  » •  ♦  *  *io 
l  n 


Schoenberg  (1938)  pointed  out  that  ♦  Z)  9^  if  m  <  n.  If  the 
distribution  function  F  of  R  is  related  to  <)>  £  4>n  as  in  (2.3) 
with  n  substituted  for  m,  then  also  <j>  £  $m,  m  <  n,  and  there 

it  it 

exists  a  distribution  function  F  (x)  of  R  being  related  to  $  as 

it 

in  (2.3)  with  F  substituted  for  F.  Cambanis,  Huang  and  Simons  pointed 
ft  d  2 

out  that  R  *  R*b,  where  b  ^  0,  b  ^  B(m/2, (n-m)/2)  and  b  is  indepen¬ 
dent  of  R.  For  convenience  we  denote  these  relationships  by  R  •*-»-  <{>  € 

“*  R‘  3  R  bm/2 ,  (n-m) /2  ~ 


X  -V  MECnXp(M;E1,...,En;$) 


with  rk(Ej)  ■  k^t  j  ■  1, 


if  and  only  if 


(2. A) 
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(V  OO  (k_) 

where  1^,...,!^  are  Independent  of  J  (j-l,...,n),  Uj^  . 

are  independent,  I  -  A! A.  is  a  factorisation  of  I  ,  j «* 1, . . . ,n,  and 
the  joint  distribution  function  F(x^,...,xn)  of  (R^,...,Rn)  is 
related  to  <p  as 

(2.5)  <Ku, . U  )  -  [  •••  [  fi  (l\U.) . (rjl  )dF(r. . O. 

J  [0,«)  >  [0,»)  kl  1  1  kn  n  n  1  n 

and  r,X  -  Y"  denotes  that  the  distribution  of  X  is  the  same  as  that  of  Y. 

The  proof  is  similar  to  one  of  Theorem  1  (Cambanis,  Huang  and  Simons 
(1981)).  The  properties  of  the  operation  are  discussed  by  Anderson 

and  Fang  (1982).  The  following  two  properties  are  important  in  this 


paper: 


(i)  Assume  that  X  *  Y  and  f. (•),  J*l,...,m,  are  Borel  functions. 


fl(X) 


fl(Y) 


f  (X) 

m  — 


f  (y) 

m  - 


(li)  Assume  that  X  and  Y  are  n*p  random  matrices,  z  is  a 
random  variable  and  is  independent  of  X  and  Y,  respectively.  If 


(2.6) 


p(z  >  0)  *  1 


then  X  *  Y  if  and  only  if  zX  -  zY. 

By  using  the  first  property  and  (2.4)  we  have 

(2.7)  (rJ,...,r£)  $  h 


•'(?(n)"i!(n)),5n(5(n)"y(n) 


5 


where  A~  is  a  generalized  inverse  of  A.  In  particular,  if  E.  ■  I  9 

-  ~  ~j  -P 

j  »  l,2,...,n,  then  (2.7)  becomes 

(2.8)  (R^,...,Rn)  ^  ^ll^(l)  11^*  *  *  *  »li*(n)  * 

2 

where  ||x^||  =  x(i)x(i)»  1  =  1 . .  or 

(2.8) '  (R^,...,Rq)  »  ^ll*(i)  II»  •  •  •  »||*(n^|] )  • 


In  this  case,  if  p(X-O)  *  0,  we  have 


(2.9) 


1?^ . 


(n)  1 


(u 


(p) 


When  X -x-  MECixp(M;E1,...,En;4»)  with  . . .  ,un)  *  exp[-(u]L+- •  *+un) /2] , 

the  corresponding  population  is  normal,  i.e.  x^^  %  Np^(j)’^j^’  ^  . n» 

2  2 

and  x(i) » • • • »x(n)  are  independent.  Now  R^,...,Rn  are  Independently  dis¬ 


tributed  according  to  Xv  »  respectively.  If  the  joint  distribution 

X1  Si 

2  2 

of  R£,...,Rq  is  a  multivariate  chi-square  distribution  (or  generalized 
Rayleigh  distribution) (cf.  Johnson  and  Kotz  (1972),  p.  220),  then 

;(j)"V^)-y'3,1’-’n’  but  5(1)  *  *  *  *  »5(n)  are  dependent.  By 
this  method,  we  could  generalize  the  theory  of  the  distribution  and  esti¬ 
mation  in  the  multiple  normal  population  to  the  dependent  case. 

d  (V 

Suppose  X  ^  MECnxp(M;E1,...,En»<J)).  From  Lemma  1  X^  *  |J(j)+Rj5j5j  ^  * 

*  * 

i.e.  ^  ECp(ii(j)  ,1^ ,  <j>  ),  where  <f>  €  R^ .  It  is  easy  to  see 

that  <|>*(u)  ■  <J>(0, . . .  ,0,u,0 . 0)  where  u  is  in  the  j-th  position. 


What  is  the  marginal  distribution  of  x. ?  From  (2.4)  we  have 

~J 
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(2.10) 


5)  '  Sj  + 


* .<»>  ,(v 

n~j  ~n 


.<« 


where  a^  is  the  j-th  column  of  A^,  j«l,...,p;  £  -  1 . n.  Let 

h  "  <°«)>  *»*  h  “  <?£i . Uk*'  *  V°’V’  then 


(£),  (V  d 

«j  «£ 


as 


*  aj  j  ^  *  and  /U  ?«.  II2  ^  B  [§••  ~~~2  »  (cf  •  Anderson  and 


Fang  (1982)).  Hence 


(2.11) 


w 


~J 


Hj  + 


o(1)  . 

Jj  *1 


I R  CJ?  Z  I 
l  n(  jj  nj 


k  -1 


2  i 

where  z^, • . . ,zn  are  independent ,  *J  B<2-  2  >•  *1 . *„ 

are  independent  of  z.  , . . . ,  z  and  (R, , . . .  ,R  )  'V*  F(x. , . . .  ,x  )  . 

in  in  in 

We  believe  some  further  results  similar  to  the  elliptically  contoured 
distribution  could  be  obtained  by  the  method  used  by  Cambanis,  Huang  and 
Simons  (1981).  In  this  paper  we  mainly  pay  attention  to  a  specific  sub¬ 
class  of  multivariate  elliptically  contoured  distributions  which  will  be 
defined  in  the  next  section. 
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3.  The  Case  in  Which  the  Characteristic  Function  is  Composed  of 
Addition  of  Arguments. 

From  now  on  we  assume 


(3.1) 


<Kun,..  .,U  )  *  <Ku.  +  •••+u  ), 
J.  n  l  n 


and  we  still  denote  the  multivariate  elliptically  contoured  distribution 
by  MECnXp(M»£i»-**»£n;<J>). 

Theorem  1.  X  'v  MEC  (M:I. .....E  ;d>)  if  and  only  if  x  'u  EC 

—  nxp  ~  „n  T  „  np 

(y»V,4>)  with 


(3.2) 


?1  0  ...  o  1 

0  I2  ...  0 


•  •  • 

•  •  * 

0  0  ...  E 

~  ~n 


Proof .  If  X  'V*  MEC  (M;E. , . . .  fE  ;<j>),  the  characteristic  function 
- — — —  .  nxp  —  —  j.  ~n 

of  X  is  (cf.  (2.1)  and  (3.1).) 


exp 


l1  “ 


exp(it'y)<f>(t'Vt)  , 


i.e.  x^EC  (y,V,$).  The  "if"  part  is  obvious  and  the  theorem  follows. 
~  np 


Theorem  1  shows  us  that  for  any  <p  €  ^  we  can  construct  a 

♦Ok  »•••»»»  )  by  (3.1)  such  that  *  €  $.  t.  with  k-  +  ‘**  +  k  ■  k. 

1  n  K.).ti|K^  l  n 
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Corollary.  X  ^  MEC  (M;Zt...,Z;$)  if  and  only  if  x  ^  EC  (y,V,4>) 

*'*  »  P  -  ^  ‘w  ^P  ** 


with  V  -  I  ®E. 

«*  *»n  **» 


Definition  3.1.  If  the  random  vector  (z . z  )'  satisfies 

■  l  m 


(3.3) 


d  2 

•  •  •  »*u|)  *  ^  (d^»  •  •  •  ,d^)  , 


.m-l 


where  (d1# . . .  .d^)  ^  . VrV*  V1"  Ei-l  di»  R  >  °’  R^F(x) 

and  R  is  independent  of  d. ,...,d  . ,  then  we  write  (z. ,...,z  )  ^ 

i  m-x  x  n 

G  (a.. , . .  •  ,a  ,,a  ;<{>)  and  (z.  ,...,z  .)  ^  G  (a.,...,a  .  ;a  ;<J>),  where 

m  i  m-i  m  i  m-i  m  1  m— 1  m 

<j)  €  $  **■  the  distribution  function  F(x)  of  R  and  n  ■  2  Em  a. . 

n  1  i 


Anderson  and  Fang  (1982)  point  out  that  if  P(R-O)  =  0  the  density 

of  zn i • • • • z  «  is 
1  m-1 


(3.4) 


,  „  m-1  a  -1 

f-(n-2)(r2_  j  >  »  dF(r)  > 

1  1 


for  z. , . . . ,z  .  >  0  . 
x  m-i 


Further  if  x  ^  ECn(y,E,<J>)  has  a  density,  which  must  have  the  form 


(3.5) 


|  E  |  ~!\(  (x-y)  *E_1(x-y) ) 


for  a  suitable  function  g(*)>  then  the  density  of  R  related  to  <t> 
is 


(3.6) 


.  2irn/2  n-1  ,  2. 
f(r)  - - —  r  g(r  )  . 

T(|) 
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In  this  case  the  density  of  G  (a,,..., a  ,  ,a  ; <h)  has  the  simpler  form 

in  x  jj-i  m 


(3.7) 


n/2  m  a  -1  m 

n  r(a,>  1*1  1 

1  1 


From  Theorem  1,  if  X  'V  MEC  (MjZ, ,...,!  ;<j>)  and  rkZ  -  k, ,  j  «  1, . . . ,n, 

~  nxp  ~  ~  i  -n  ^-J  J 

then  x  'v*  ECn^(y ,V,<j>)  where  V  is  defined  by  (3.2).  Hence,  x  has  the 
stochastic  representation  (cf.  Cambanis,  Huang  and  Simons  (1981)) 


(3.8) 


U  +  KB'u(k)  ,  k  =  lK, 


where  B'B  *  V  and  B'  is  an  np *  k  matrix.  On  the  other  hand,  there 
is  the  stochastic  representation  (2.4)  for  X;  what  is  the  relationship 
between  R  and  R^,...,Rn,  and  between  B  and  (i«l, . . . ,n) ? 

fit  It  \a  > 

Theorem  2.  (i)  (R^ . R*)  *  Gn( ~ . -S— ,  y-;4>  ,  (ii) 


(3.9) 


6i  9  •  •  •  o 

b  **  9  62  •••  9 


0  0  ... 


Proof .  Comparing  (2.4)  and  (3.8),  the  formula  (3.9)  follows  and  we 


r.  <Vi 
Vi 


,(H  i  : 


R  uM 

n  n  ' 
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Let  ?kxi  *  VS'Jk5  y  “  (y(1)  . y(n^  >' »  Where  y^,...,y^n* 

have  1t^ , . . . , elements  of  y,  respectively.  Therefore 


|y(1) 

I|y<1)ll1 

’  lly(l)ll 

« 

(kl) 

Ru(k)  4  R 

E  y  II 

• 

• 

• 

-  R 

lly(1> 

II  II  y  II 

• 

• 

d 

'hip- 

Hi 

(n) 

y 

l|y<n)ll 

l|y(o)ll 

<V 

|y|| 

. 

' 

l|y<n) 

. 

II  (|y  II 

* 

RT7IT 

u 

4 

because  ^  «  2^/||y^||,  j  *  lf...,n,  and  y^/||y^||  is  independent 

of  ||y^^||,  and  therefore  y^/||y^'j|  is  Independent  of  ||y^||/|y||.  Now 


d 


1yfl)  Ip 
.  IIP2 


» 


the  theorem  follows.  Q.E.D. 


Corollary  1.  If  ^  “  l2  "  *  *  *  "  Zn  “  l  *  ~'6»  rk(P  “  £»  then 
(^.....R*)  *  Gn(£/2,...,*/2,fV2;4>)  and 


B  - 


A  0  ...  0 

0  A  ...  0 

•  •  • 

•  •  • 

0  0  ...  A 


Definition  3.2.  If  X 't  MEC  ..  (M;E, , . . .  ,E_  ;<J>)  ,  E,  -• 
"""  11  ^  nxp  -  ~l  ~n 

end  M  ■  e  ®p',  we  write  X 't  LEC  „  (y,E,4»). 

-  ~n  ~  —  n*p  ~  ~ 


E  -  E 
~n 
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Corollary  2.  Assume  X  'v*  LEC^Cy, £,<{>)  and  rk(E)  ■  l,  then 
(3.10)  X  -  en®  y’  +  RUA  , 

where  A'A  ■  £,  A:  8,*p,  U:  nx£,  vec  U  -  ,  R  **  $  e  and  R 

Is  Independent  of  U. 

Proof.  From  Theorem  1  and  Theorem  2  we  have  x  ^  EC  (e  ®y,I  ®£,<J>). 
—  -  —  -  ~~np  ~n  ~  ~n  ~  ~ 

Let  U«  (u(1j . u^n))  '  ,  then 


fA'  ...  0  I 


x  2  e  ®  y  +  R 

~  — H  -s. 


0  ...  A’j 


(nJt) 


e  ®y  +  R 

~n  ~ 


'i'ia)' 


'-(n)  J 


.  (nf,) 

because  u 


vec  U 


d 


vec  U' . 


Thus 


x  2  e  ®ii'  +  R 


~(n)~- 


E  <2>fl'  +  RUA  .  Q.E.D. 

n  ->» 


In  order  to  obtain  the  marginal  distribution  of  X,  we  need  the 
following  Lemma  which  is  from  Cambanis,  Huang  and  Simons  (  1981),  and 
Anderson  and  Fang  (1982) . 
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Leona  2.  Assume  y  EC  (y,I  ,d>),  4>  €  $  corresponds  to  R  and  C 
■  -•■■■  —  ~  n  ~  ~n  n 

is  an  n  x  p  matrix  with  rk(C)  -  p  <  n,  then  x  -  C'y  ^  • 

where  <0  €  *p  ~R*  -  »p/2><n-p)/2- 


Assume  X  ^  LEC  (y,E,<|>),  rk(E)  ■  i  and  <f>  €  •*+  R,  from  Theorem 

~  n  p  ~  ~  —  n* 

1  and  Lenma  2  we  immediately  obtain  the  marginal  distributions  of  rows 
and  columns  of  X: 


(1)  ?(j)  *  where  *  €  ~  R*  -  RbA/2,£(n-l)/2- 

(2)  ?j  %  BCn(pjen,ojjIn,*),  where  $  *  *n  ~  **  =  Rbtl/2,n(A-l)/2’ 
is  the  j-th  component  of  y  and  E  =  (cr^). 


The  following  corollaries  concern  the  distributions  of  linear 
functions  of  X. 


Corollary  1.  Assume  X^  LECnXp(y ,E ,$) ,  rk(E)  *  l,  ♦  €  *nl  R» 
and  B  is  a  p  *  q  matrix  with  rk(B)  =  min(p,q),  then 

(1)  XB  %  LECnXq(B'y,B'EB,4))  with  <i>  6  R  if  q  1 

(2)  XB  *  LECnXq(B’y,B'EB,<J»)  with  4>  €  *nq  ~  R*  *  Rbnq/2,n(Jl^)/2 
if  q  <  JL 


Proof.  As  X  *  LEC__(y,Z,<J>)  with  rk(E)  -  l,  then 

■  <v  nxp  v  ^  ~ 

X  -  e„<g>  y’  +  RUA  , 

<v  vll  <v 


where  A' A  -  E.  Thus 
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or 


-  eQ  0  (y’B)  +  RU(AB)  , 

B'A' 

vec(XB) '  -  @  (B'y)  +  R 

B’A' 

The  corollary  follows  from  Leona  2  and  Corollary  2  of  Theorem  2. 

Corollary  2.  Assume  X  LECnXp(p,E,^) ,  rk(E)  ■!,<!>€ 
and  B  is  a  q  *  n  matrix  with  rkB  ■  q  <  n.  Then 


(3.11)  vec(BX)  'V.  EC  (y®(Be  ) ,  E0(BB,),<f>)  , 

qp  ~  ~~n 

"h“*  *  e  V  ~R*  ’  %m.(n-qn/2- 

Proof.  From  the  assumption 

X  -  en®  y'  +  RUA  , 


Thus 


vec(BX)  ■  y  ® (Be  )  +  R 

~~  —  — .n 


=  y©  (Be  )  +  RC’u(n£)  (say)  . 

"  ~~n 

Using  Lemma  2  vec(BX)  ^  EC  (y©(Be  ),C'C,d>)  with  d>  £  ^  „  ++ 

~~  qp  ~  ~~n  -  -  qx, 

R*  *  Rbq£/2  (n-q)Z/2‘  The  coroolary  follows  from  C'C  *  E  ©  (BB ')  .  Q.E.D 

Remark.  The  above  two  corollaries  show  us  that  the  distribution  of 

XB  still  belongs  to  the  class  of  the  multivariate  elliptically  contoured 

distributions,  but  the  distribution  of  BX  in  general  does  not  belong  to 

2 

this  class  unless  BB'  =  c  I  where  c  is  a  constant. 

~~  ~q 

Theorem  3.  Suppose  that  X  'v  LECnXp(y  ,E,<|>)  with  E  >  0,  then 

(1)  Rz  =  tr  E  G  where  R  ■♦+■  $  €  $  and 

~  ~  np 

n 

(3.12)  G*  J^(X(i)-p)  (X(i)-y)  '  *  (X-en©y')  '(X- en©  y’)  ; 

(2)  The  density  of  X  has  the  form  of 

(3.13)  |Ep/2  g(tr  E_1G) 
if  it  exists. 
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Proof.  As  X  'v*  LEC _ (jj,E,  <#>)  ,  then  x  “  vec  X*  ^  EC__(e_  ®y,I_  <S>  E ,  <f>) . 

1,1  ~  n*p  ^  v  ^  —  np  'vii  %  vii  -v 

From  Corollary  1  of  Cambanls,  Huang  and  Simons  (1981) 

R2  2  (x-en  ©  p)  ’  (ln  <2>  p_1(x-en  ®  p) 

-  (x-en®  y)  '(In®  E_1)(x-£n®  P 

■  tr  E  . 

If  X  has  a  density,  so  does  x  and  the  density  of  x  has  the  form  of 
(cf.  (3.5)) 


|l  ®  El-1*  g((x-e_®  U)’(I  ®  E)-1(x-e  ®  y)) 
‘  ~n  ~  ~  ~n  ~  ~n  ~  —  ~n  ~ 

”  \l\~n/2  g(tr  E_1G).  Q.E.D. 


Corollary  1.  Under  the  supposition  of  Theorem  3,  if  X  has  the 
density  (3.13),  then  R  has  the  density 


(3.14) 


2TT>*np  rnp-l 
T  (^  np) 


g(r2) 


The  proof  uses  (3.6). 


Corollary  2.  Assume  that  x^  ^  Np(jj,E)  with  E  >  0,  j  n  and 

-1  2 

X(l) *  *  *  * *x(n)  ar*  ^dependent,  then  tr  E  G  ^  X  np* 
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2  2 

Proof.  Take  $  -  exp(-t/2)  in  Theorem  3.  As  we  know  R  np> 
the  corollary  follows. 


4.  The  Distributions  of  Correlation  Coefficients,  the  Multiple 
Correlation  Coefficients,  and  T^ . 

Throughout  the  rest  of  this  paper  we  assume  X  ^  LEC  „  (£i,E,<£)  with 

—  nxp  *v  *• 

E  >  0  and  €  ♦  •**  the  distribution  function  F  of  R. 

-  np 

According  to  the  common  definition  of  the  correlation  coefficient  r^ , 

(;r*  A.1 ' 

where  x  *  —  e'x. ,  i  -  l,...,p.  Let  D  *  I  e  e';  then 
i  n  ~n~i  ~n  n  ~n~n 


(4.1) 


x  Dx 


^  /x’Dx.  •  x'Dx 
Y  ~  i.  -.j 


At  first  we  consider  x/4x  'v*  N  (y,E),  j  *  l,...,n,  and  suppose 

~vj/  P  ~  ~ 

x,,.,...,!.  .  are  Independent.  There  exist  y,..  ^  N  (0,1  ),  j  =  1 v • . 
~vjJ  ~vn;  ~Vj7  P  ~  ~P 


^(l)**"’2(n)  Independent  such  that 


(4.2) 


*ii)  ’  ~  +  -’y(j)  ’  * 


where  E  ■  A' A.  Let 


2(1) 

H 

(4.3) 

Y  - 

• 

• 

,  A  “  (a.,..., a)  and  J  • 

• 

• 

l\nV 
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it  is  easy  to  see  that 


(4.4)  +  I?i  *  1  "  1 . .  • 

Hence  we  have 


(4.5) 


ajY'DYa^ 


j*  m  ^  -i 

iJ  (a'Y'DYa.^U.'Y'DYa.)5* 

~  -J  —  — —J 


Secondly,  we  come  back  to  the  case  of  X  ^  LBCnX^(p,H,<p) .  From  (3.10) 


~j  "  yj~n  +  R~~J  ’ 


where  U:  n  *  p  and  vec  U  =  u^n*^  .  Hence 


x'Dx,  -  R2a!U’DUa,  -  R2a!Y'DYa4/tr  Y'Y,  V  i,j, 

"*»  l~*wj  •»  1**  v<v<wj  «•  1*»  - - -  1  *v 


because  U  •  Y/(tr  Y'Y)1*,  now  (4.1)  becomes 


.  a'Y'DYa. 

ij  (a'Y'DYaJ^a.'Y'DYa.)^ 


which  is  equivalent  to  the  normal  case  by  the  properties  (i)  and  (ii) 
"d" 


of  the  operation 


Theorem  4.  Suppose  that  X  'v  LEC  v  (y, £,<{>)  with  £  >  0,  then  the 

~~  «*•  n*p  *«*  •» 

joint  distribution  of  r^,  1  -  l,...,j-l;  j  -  2,...,p  is  the  same  as  the 
normal  case  where  x^  ^  N^d),!),  j  -  l,...,n,  and  x(i) » • • • >x(n)  are 
independent . 


Corollary  1.  Suppose  that  X  *v*  LEC  (y»l  ,4»)  and  R  -  (r,.), 

*  n*p  —  ~p  —  xj 

then  the  joint  density  of  r ^ ,  i  <  j  is 


(4.6) 


[r(-2“)]  |R|W»-P-1) 

rp<7">  ' 


where  V  (?)  *  ^  n  r (—  S*— )  and  m  *  n-1.  In  particular,  the 


density  of  r^  is 


(4.7) 


r(2n)  (1  2  A(m-3) 

/F  r(|(m-l))  (  ^ 


(cf.  Anderson  (1958),  Section  7.6). 


Corollary  2.  Suppose  that  X  ^  LEC  (il,E ,4>>  and  p  *  a  //j  c 

~  nxp  ~  ~  I]  ij  ii  Jj 

j*  0,  then  the  density  of  is 


(4.8) 


2m'2(l-p2.)mA,J(l-r2  )(m"3)/2  »  (2p,.r 


(m-2) !  tt 

(cf.  Anderson  (1958),  Section  4.2). 


r2(|(B4a))  . 


Now  we  consider  the  distribution  of  the  multiple  correlation  coefficient, 
for  instance,  the  multiple  correlation  coefficient  between  the  first 
variable  and  the  rest  of  the  variables.  Denote  (cf.  (3.12)). 


G  - 

811 

CM 

H 

and  E  “ 

'cn 

Su' 

Su 

-”J 

1 22 

where  g..  and  o..  are  the  first  diagonal  elements  of  G  and  I, 
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respectively.  According  to  the  definition  of  the  sample  multiple  correla¬ 
tion  coefficient 


r: 


2  5l2?22521 


g 


11 


It  is  easy  to  see  that  G  “  X’DX,  •  xiDxl  •  ~12  -  921  "  DX^2),  and 

G22  “  X'(2)DX(2)  where  X(2)  "  (x2»  •  •  •  *xp^  *  By  the  same  meth°d  as  1x1 
the  case  of  the  correlation  coefficients  we  have 


ffia a 

a^Y’DYa^ 


where  A...  ■  (a. )  for  the  normal  case  and  the  same  expression  for 

~\2)  ~2  -p 

2 

R.  in  the  case  of  X  ^  LEC  w  (£,£»$>)• 

w  n^p  -v  ~ 


Theorem  5.  Suppose  X  LEC  .._(C,E,4>)  with  E  >  0. 
“  n*p  —  ~  ~ 


(1)  If 


-2  _  ~12~22~21  „ 

R  =  - - - 0  , 

11 


then  [R*/(l-R^)](n-p)/(p-l)  ^  F(p-l,n-p) 


(2)  If  R2  >  0,  then  the  density  of  R2  is 

(l-R^)J*(n"P‘2)(l-R2)Js(n_1)  »  (R2)C,(R2)5s<P"1)'fCt‘:L  r2(-j(n-l)+a) 

T(|(n-p))r(|(n-l))  a-0  at  r(|(p-l)4a) 

(cf.  Anderson  (1958)  Section  4. A)). 
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Suppose  X  LEC  (y,E,$)  with  E  >  0  we  want  to  test 

*v  n*p  *%»  ***  •• 

H^:  {j  ■  Pq  and  H^:  y  4  y^  . 

Without  loss  of  generality  we  can  assume  y^  -  0.  It  is  well-known  that 

2 

we  can  use  Hotelling's  T  -test  for  testing  the  null  hypothesis  under  the 

2 

normal  distribution.  What  is  the  distribution  of  T  in  the  case  of  the 

multiple  elliptically  contoured  distributions? 

2 

According  to  the  definition  of  T  we  have 

(4.9)  T2  ■  n(n-l)x'G~^x  , 

where 

<4-10>  • 

and  G  is  defined  by  (3.12).  Using  G  *  X'DX  and  (4.10)  we  have 

(4.11)  T2  -  n(n-l)(i  e£) (X'pX)'1^  X'en)  . 

In  the  normal  case  we  have  X  -  e  ®  p'  +  YA  with  A' A  ■  E  (cf. 

~  ~n  -  ~  —  -  —  - 

(4.2)-(4.4)) .  If  the  null  hypothesis  is  true  we  have 

(4.12)  t2  -  —  e,YA(A’Y'DYA)"1A,Y'e  . 

n  ~n —  ~  ~  — ~  ~  ~  ~n 

By  (3.10)  and  U  -  Y/(tr  Y'Y)*5,  when  the  null  hypothesis  is  true  we  obtain 

*>»  <*•  <N* 

the  same  expression  as  (4.12)  for  T2  In  the  case  of  X  ^  LEC  (p,E,$). 


Theorem  6.  Suppose  that  X'VLEC  .  (u.E.di)  with  E  >  0  and  T2 

"  -  nxp  '  ~  -v 

Is  defined  by  (4.9),  if  p  «  0  the  distribution  of  [T2/(n-l)][(n-p)/p] 
is  F(p,n-p) . 

5.  The  Distributions  of  the  Sample  Covariance  Matrix  and  Generalized 
Variance. 

5.1.  The  distribution  of  the  sample  covariance  matrix.  It  is  a  well- 
known  fact  that  the  distribution  of  the  sample  covariance  matrix  of  multi¬ 
variate  normal  population  is  the  Wishart  distribution.  If  the  sample 
X(i),...,x(N)  is  from  the  population  of  Np(p,E)  and 

N 

*  ■  Ji  • 

then  A  ^  e"  .  z  z'  where  n  -  N-l  and  z,,. . z,  .  are  l.i.d. 

distributed  according  to  N^(0,E).  Now  we  want  to  find  the  corresponding 

distribution  for  multivariate  elliptically  distributions.  Assume 

X  'u  LEC  (0,E,4>)  with  E  >  0  we  want  to  obtain  the  distribution  of 
n*p  «v  —  <v 

(5-1>  S  ■  *(«,;<«,  •  i'i  -  <V  ■ 

(1)  Assume  X  has  a  density.  Theorem  3  shows  us  that  the  density  of 

X  is 

(5.2)  |E|'n/2  g(tr  Z’^'X)  -  |Ern/2  g(tr  E_1W) 

as  p  -  0  and  W  ■  G.  By  Lemma  13.3.1  (p.  319,  Anderson  (1958)),  the 

density  of  W  (l.e.  the  density  of  w.,,...,w,  ,w.„ . w.  , ...,w  )  is 

~  u  ip  it  zp  pp 
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(5.3) 


7 


np/2-p(p-l)/4  /  i  \  /  <)  _ / o  ^ 

■  |wj  <n“p-1>/2  |r|"n/2  g(tr  E_1W) 


n  rc 

a-l 


n-a+1 


) 


In  the  normal  case  (5.3)  reduces  to  the  Wishart  density. 


(2)  Assume  p(X=0)  «  0.  In  this  case  it  is  possible  that  X  does 
not  have  a  density,  but  all  the  marginal  distributions  of  X  will  have 
densities  (cf.  Kelker  (1970)).  Ue  consider  the  following  interesting  and 
Important  situations: 

(A)  Let 


(5.  A) 


*k 


'll 


nl 


lk 


nk 


and  *  X^  ,  1  £  k  <  p 


In  order  to  obtain  the  density  of  X^  we  need  the  following  Lemma. 


Lemma  3.  Assume  x  ^  ECn(0,Z,4>)  with  Z  >  0  and  p(X=0)  *  0. 
Let  x  ■  w^ere  x(i)  is  an  m  *  1  vector,  m  <  n,  then 

the  density  of  x^  is 


(5.5) 


r(|) 

V^rc 


r 

'(^r>  J(xJ 


£(l)£llJ(l)> 


-(n-2) ,  2  ,  j-1  ^ (n-m) /2-1  dF(rj 

J*  r  U  X(l)  11  (1);  ' 


where  Z^  is  the  first  principal  minor  of  order  m  of  Z. 


Proof .  If  Z  ■  In,  the  density  of  x^  is 


23 


r 


(5.6) 


T(f) 


TTm/2  rc 


r^r)  ^x'  .X 


(1)~(1) 


r"(n"2)(r2-x*1)x(1))(n^l)/2-l  dF(r) 


(cf.  (21),  Cambanls,  Huang  and  Simons  (1981);  the  lower  limit  of  the 
Integral  In  that  paper  should  be  >/u,  not  u) .  When  E  j  I  the  Lemma 
follovs  by  using  the  transformation  y  *  E~x.  Q.E.D. 


4 


Using  Lemma  3  and  noting  the  structure  of  X.^,  the  density  of 
is 


(5.7) 


r  <rf)\z 


l-k/2 


2  "'tk1 


^nk/2  r(ni£zkl)  J 


r 

i . - 

/  trir1X1' 


r-(np-2)(r2_trJ;^lx^x^)n(p-k)/2-l  dp(r)  ^ 


k  *k?k 


where  Efc  is  the  first  principal  minor  of  order  k  of  E.  By  Lemma 
13.3.1  of  Anderson  (1958)  again,  the  density  of  W^  is 


(5.8) 


r(I?)  l~k!"n/2 


Trk(k-l)/4r(rL(p_k))  n  r(n=|fl) 
a*l 


«?k 


(n-k-l)/2 


-(np-2) 


(r2-trE^1Wk)n<P‘k)/2“1  dF(r) 


(B)  Let 


X 


(m) 


xn  x, 
11  . . .  Ip 


x  1  ...  x 
ml  mp 


p<m<n,W,  .»X!  NX,  . 
*  -  ~(m)  ~  (m)  ~  (m) 
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are  found  as 


In  a  similar  way  the  densities  of  X,  ,  and  W,  . 

~(m)  ~(m) 


(5.9) 


r(f)|zr/2  p 

^p/2  T(^h.) 


/trE_1X) 


r-(np-2)(r2_trE-lx, 

~  ~(m)~(m) 


(m)?(m) 


and 


(5.10) 


r(^)|ifm/2 


7rp(p-l)/4r(ln^£)  K  r(m=|tl) 
a*l 


lw  [  (m-p-l)/2 
1  ~(m) 


-(np-2) (r2_trE-lw  )(n-m)p/2-l 
~  ~  W 


respectively. 

Further,  partition  X  into  k+1  parts,  i.e. 


'*L 

h 

ifi 


where  X^»X2»  •  •  •  are  x  P»  n2  x  p’,***nk+l  x  p  matrices, 

k+1 

tively  and  p  _<  n.^  <  n,  i  *  l,...,k,  1,  E^  n^  ■  n.  Let 

tt,..  *  X'X.,  i  *  l,...,k,  then  the  joint  density  of  W.  . . W.  . 


^ (n-m)p/2-l 


dF(r)  , 


respec¬ 


ts 


(5.11) 


-(n-n  )/2 

r(^)|z|  k+1 


kp(0-iy4  nk+ip  k  p  ni*^+1  i-i  ~(J) 

tt^p^p  y  rc— ^-)  n  it  rc-J  .  >  J  1 

j«l  0-1 


k  (n  -p-l)/2 

n  w,.J  j 


f“(np-2)  fr2  _  ]^l/MdF( 

k  -1  I*  4,1  ~  ~(j>J 

(A tr^  3 


The  mechod  to  obtain  (5.11)  Is  the  sane  as  to  obtain  (5.10),  but  here  we 
need  to  use  Lemma  13.3.1  of  Anderson  (1958)  k  times.  When  p  -  1, 
the  density  (5.11)  reduces  to  Gk+1(n1/2 , . . . ,nk/2;nfc+1/2 ;$) .  We  denote 
the  density  (5.11)  by  MGp  k+1(Z;n1/2, . . .  ,nk/2;nk+1/2;«J>)  . 

Further,  if  X  has  a  density  (3.13),  we  rewrite  it  as 


fk+1  , 

m-  »1  It  trrt(j) 


We  use  Lemma  13.3.1  of  Anderson  (1958)  k+1  times  to  obtain  the  density 

°f  ?(!)»* •••?(k+i)  as 


n  n  /  0  _  .  _ n  /  O 


(n  -p-l)/2 


n.-a+l 


r 

f 


2iTnp/2-kp(p-l)/4|I.|"(n"nk+l)/2  k 

-  n 

j-i 


r< 


nk+ip  k  p  IN"a+1 

\  TT  TT  n  /  yi 


-)  n  n  r(- 

J-l  a-1 


) 


f 


!(j)' 


(nj-p-D/2 


nk+ip-1  r 
*  * 


r"+  l  trI_1VJ 

J-l  '  “ 


(j) 


dr. 


(C)  Consider  the  marginal  density  of  p-l»wip» *  *  * *wi  p_i» 

w2p . w2,P-i*  *  *  *  ,wp-2,p-i,wip . Vi.p  (i'e‘  the  aar^inal  density 

of  W  except  wpp);  f°r  simplicity  we  will  say  the  density  of  W  . 


(1)  Let  E  *  U'U  where  U:  nxp  and  vec  U  *  u 
defined  by  (4.3),  then 


(np) 


Let  Y  be 


(5.13) 


E  =  Y'Y/trY'Y  -  (e^)  . 


We  want  to  find  the  density  of  E  (i.e.,  the  density  of  E  except  e  >. 

**”  PP 

However,  the  correlation  coefficient  r^  is 


A 1 


,/^J  e  /£  f 
ii  jj  ii  jj 


where  (f^)  =  F  =  Y'Y.  As  we  know,  {r^  ,  i*  1, . . .  ,J-l,j-2 . p}  are 

independent  of  (JlyJI,  i-l,...,p}  and  (ilyx I|2 » •  •  •  *!l yp_i ll*"5 /tr^ ’Y  % 
Dpt|.  •  •  •  ,y; j)  ,  hence  (e^ , . . . ,  ep  l  ^  p  l)  %  Dp(y, . . .  ,^;^)  and  the  density 


of  e, 


eH>  *  •  -  *e. 


,c  .  -  /e  "  7  is 

’  P-1  P-1 » P-1 


>P-lr 


(5.14) 


V  e-1 


i-1 


1  '  ?l  \ 

i-1  1 


f  ■ 1 


Pr1  2 

if  I  e2.  <  1  . 
1  1 
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From  Corollary  1  of  Theorem  4,  the  joint  density  of  e^,...,ep  i»ri2****» 
rlp’  r23’  *  rp-l,P  18 

2pr® 

(5.15)  - 2- - 

^(p-D/4  „  r(nriti) 

i-1  l 

i=l 


l«|«T«  V  ,i.fA 

i-1  1  1 


prl  ,  j"1 


P-!  2 

if  l  e.  <  1 


Consider  the  following  transformation 


2 

ei 

eiejrij 


i  -  1 , . . . , p-1 
i  -  1, • . • » j-1;  j  - 1» • • • pp> 


with  e 


P“1  2  k 

-  (1  -  l  e. )  .  It  is  easy  to  see  that 
1  J 


the  matrix  of 


3(e 


11’ 


Vl,p-l’c12’  •  •/  ’Xp-l’e23”  •  •  ’e2lP-l . ep-_2Ji>-l!elj> 


» •  •  • » e_  i  _) 


9(*1’  *  •  *  •' Vl’r12*  *  *  *  ,rl,p-l*r23*  *  *  •  ,r2 ,p-l*  *  ■  *  ,rp-2,p-l’rlp . rp-l.P> 


is  lower  triangular  with  diagonal  elements  2e^,.  . .ie^j^e^ej, • . • »eiep_i*c2e3* 
. ...e.e  . .  „e  , ,e, e  , e  ,  thus  the  Jacobian  of  the  transfor- 

'  2  p-1’  ’  p-2  p-1’  1  p  p-1  p’ 

mat ion  is 


(5.16) 


2p-l  e(p-l) 
P 


p-1 

n 

i-i 


(p-l)/2 


p-1 


i-1 


Combining  (5.15)  and  (5.16),  the  density  of 


* 


E 


18 


28 


(5.17) 


r£E) 


*p(p-1)/4  n 
1-1 


Wn-p-l)|x.Pj1eu|'i(n'P'1)  ^  ^(n-P-D 


i-1 


11 


r(^) 


^(p-D/4  P  r(nZ|±l) 
1-1 


I  Js(n-p-l) 


with  e  -  1  -  E?  1  e  . . 

pp  1  11 


(2)  Let  E  be  pxp  a  positive  definite  matrix  and  A  be  an  upper 

triangular  matrix  such  that  A'A  -  E.  Let  V  -  A'EA  -  (v^).  We  want 

* 

to  obtain  the  density  of  V  (i.e.  the  density  of  V  except  vpp) • 
Partitioning  V,  A  and  G  as  follows 


rVll  *1) 


l~U>  VPP 


(A*  0  ' 


til 


1^(1)  VJ 


fn  *(i) 


15(1)  epp 


~11  ~(1) 


0  a 


PP  J 


we  have 


(  v 

-  A'  E  A  . 

~11 

~11~11~11 

~<1) 

”  tiihita) 4 

V 

-  a^  e  +  2a 

PP 

PP  PP 

As  e  and  v  are  not  independent  variables,  the  Jacobian  of  trans- 
PP  PP 

formation  is 


29 


J(V  -  E)  -  J(VU  -  Eu)  •  J(v(1)  -  e(1)) 


J<»ii  *  in*  •  l*ii*ii 


I18"  ■  IJnlp/2 . 


where  E^  is  the  first  principal  minor  of  order  (p-1)  of  E,  and 

J(r<i)  *  la)’  ’  liiil^1  ■  lJul%  C  • 


Thus 


(5.18)  J(V  -  E)  *  |E11|,S(P+1)  aP_1  -  |E11[!s(p+1) 
~  ~  1 —11  *  pp  1 ~X1 1 


|E|  Wp-1) 


'111' 


|Enl|E 


Noting  |v|  -  I  A' I |E| |a|  -  |e| |E| ,  the  density  of  V  is 


(5.19) 


r(^) 


„p(p-i>/4  j  r(stai) 
i-1  Z 


— *srt+i  |  £  i  i  | —1  |y|,s(n'p_1) 


As  v  is  not  an  independent  variable,  we  need  to  find  the  relationship 
between  it  and  the  other  variables.  Since 


U'U  -  E  -  A'^A-1  , 


we  have 


1  -  trU'U  -  tr(A,-1VA~1)  -  trVE*1  . 


Denote  E  1  -  (ci:*)  and  V  -  (v^);  then 


^(p-l) 
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(5.20)  v 


PP 


-<°pVl(l-  ?  Pf  a\.-?  A.  1  :  („,)  . 

I  i-1  j-1  «  i-1  ip 


(3)  Assume  X  •u  LEC  „  (0,Z,<Ji)  with  Z  >  0  and  P(X-O)  «  0.  From 

^  n*p  ~  ~  ^ 

(3.10) 

X  -  RUA  , 

where  U:  nxp,  vec  U  ■  u^np\  A  is  a  pxp  upper  triangular  matrix  and 
A' A  -  Z.  Then 


(5.21)  W  =  X’X  »  R2A'U'UA  *  R2A'EA  -  R2V  . 

We  have  obtained  the  density  (5.19)  of  and  it  has  the  form  of 

where  c  is  a  constant.  Thus  the  cdf  of  for  W  positive  definite 

is 

2  2 
p{R  vjj  -  wjj’  1 . .  <w  ,  i<  j} 

w  w 

P{v  —  ~ ^2 ’  j“i»****P_1»  £  ~^2»  1  <  j  >  V  positive  def inite}dF(r). 

JJ  r  J  r 

When  the  probability  is  written  as  the  definite  integral  of  the  density 
of  VA  differentiation  yields  as  the  density  of  WA 

r-p(p+l)+2 | H | S(n-p-l)  dp(r)  # 


(5.22) 


r 

Jo 
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where  H  -  (h^),  h±i  -  w^/r  ,  i  -  1 . p-1;  h^  -  w^/r  »  1  <  J  and 


h  -  (opp)_1(l-w*/r2)  *  (oPPr2)  1(r2-w*)  , 
PP 


where  w  *  E P  ,  Ep  ,  a^w. .  +  Ep  ?■  a*Pw  .  As  h  must  be  positive, 
i-1  j-1  ij  i-1  ip  pp 

the  density  of  W*  is 


(5.23)  c  r  r 

J  /w* 


-p(p+l)+2  -p(n-p-l)  |  ,5s(n-p-l) 


dF(r) 


c  f  r-<pn-2)|w|l5(n-p-1)  dF(r) 


r(&) 

J  [.(Szi+i) 

i-1 


|r|Jsn+1|E„rl  [  r-(lm-2>  |w|*s<"-p-1> 

^  "'ll  ^ 

✓w* 


where  W  -  (w^)  with 


(5.24)  w  - 
PP 


(app)-1(t2-w*)  -  topp)-1fr2-  5  Pf  <>«*.. -"f  °lp»J 

i-1  j-1  i-1  Pj 


Now  we  can  summarize  the  above  results  as  follows: 


Theorem  7.  Assume  that  X  'v  LEC_„_(0,E,$)  with  E  >  0  and  n  >  p; 

nxp  V  ^ 

X  is  partitioned  into  k+1  parts  , . . . , which  are  n^xp,... »nic+i x 

k+1 

matrices,  respectively,  p  £  n^,  i  -  l,...,k+l  and  E^  n^  -  n;  W  -  X’X, 
W/.v  ■  XjX.  and  W.  is  defined  by  (5.4). 


(1)  If  X  has  a  density  (5.2),  then  the  density  of  W  is  (5.3) 
and  the  joint  density  of  w(i) » • • * »w(k+l)  i8  (**^2)* 

(2)  If  P(X-O)  *  0,  then  there  exists  the  marginal  density  of  any 

(proper)  subset  of  elements  of  W.  In  particular,  the  density  of  W^  is 
(5.8),  the  density  of  W*  (i.e.  the  density  of  W  except  w^p)  is  (5.23) 
and  the  joint  density  of  » •  •  •  »^(k)  *s  (5.11). 


Corollary  1.  Assume  X  ^  LECnXp(0,E,4>)  with  E  >  0  and  n  >  p, 
and  X  has  a  density  (5.2).  Let  4*  *  E  ^  and  V  ■  W  \  then  the 
density  of  V  is 


(5.25) 


T!_li — — IL-  |^Jn/2|v|-(n+P+1)/2  g(trW“1)  . 


n 

a=l  L 


Proof .  From  (5.3)  and  the  fact  of  the  Jacobian  of  the  transformation 
W  -  V"1  is  | v]  the  corollary  follows.  Q.E.D. 

When  <J>  »  exp(-t/2),  i.e.  X  is  from  a  normal  population,  (5.25) 
reduces  to  the  inverted  Wishart  distribution. 

Corollary  2.  Assume  the  condition  of  Corollary  1  holds  and  E  ■  1^, 
then  the  density  of  the  characteristic  roots  X^  >  •••  >  X^  of  W  is 


(5.26)  ir 


(n+l)p/2 


-p  x(»-P-«/2- 

mi  r(£f%r<*=fii». 


**1  (xi'V  s(f  xi) 
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Proof.  From  (S.3)  and  Theorem  13.3.1  of  Anderson  (1956),  the  density 

of  X, § .  •  •  ,  A  is 
1  P 


np/2-p(p-l)/4 


P 

n 

0-1 


r<£=2±i) 


p 

n 

i-i 


. (n-p-l) /2 
Ai 


^pCp+l)/* 

n 

a-l 


n  a  -x.) 

i<  j  1  J 


and  the  corollary  follows.  Q.B.D. 


5.2.  The  distribution  of  the  generalized  variance. 


Theorem  8. 
£  >  0,  p(X-O)  - 
distributed  as 


Assume  that  X^LEC  „  (0,I,4>)  with  $  €  $  *-*■  R, 

~  n*p  ~  ~  np 

0  and  n  >  p.  Let  y^,...,yp+^  independently 
VVl . Wl'Vp-l)/l'  respectively.  Then 


(5.27) 


Proof .  From  the  assumption  X  ■  RUA  **  KYA//trY*Y,  where  Y  is 
defined  by  (4.3),  then 


(5.28)  |W|  -  |X’X|  -  R2p|A’Y’YA|/(tr  Y'Y)P  -  R2p|£| |Y’Y|/(tr  Y'Y)P  . 

Let  T  ■  (t..)  be  the  lower  triangular  matrix  such  that  TT'  ■  Y’Y. 

—  ij  ~~  ~  ~ 

It  can  be  shown  that 

(a)  t^  K(0,1)  for  any  i  >  j; 
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(b)  tii,v'xn-i+i  (i-e-  tii<v,4-i+i)» *" 1 . . 


(c)  {ty,  J  <,  i)  are  Independent. 


(cf.  Johnson  and  Kotz  (1972)).  Thus 


|W|  -  R2p i Z | |TT* | /(tr  TT')P 


As  T  is  a  lower  triangular  we  have 


P  o  a  p 

|tt'|  -  n  t.  -  n  y. 

l-i  1  i-i 


and 


tr 


TT*  “  f  ‘ii  +  I  t11  “  ^  yl  +  yp+l 
lil  1  >  j  1-1  1  P*1 


P+1 

L  V 


1-1 


Q.E.D. 


6.  A  Multivariate  Analogue  to  Cochran's  Theorems. 

In  Section  5  we  obtained  the  distributions  of  and  denoted 

it  by  MG  _(£:  <J>)  .  Let  D  be  an  n  x  n  symmetric  matrix.  We 

p,2  ~  l  2 

want  to  know  a  necessary  and  sufficient  condition  for  X'DX  ^ 

*P#2(S*  f;  ^T-’  ♦>* 

Theorem  9.  Assume  that  X^LEC _ (0,£,4>)  with  £  >  0  and  P(X-O)  -  0, 

■ "  —  n*p  -n.  —  *-  - 

and  D  is  an  n  *  n  symmetric  matrix,  then  X'DX  ^  MG^  2^»  §■*  4>)  If 

2 

and  only  if  D  -  D  and  rk(D)  -  a. 

Proof.  Assume  X'DX  ^  MG  Ali  ^5  $) .  Since  X  -  RUA,  we  have 

-  ~  —  ~  £  l  ~ 

(cf.  Theorem  7) 

R2A’U'DUA  ^  X'DX  0  rVu'^A  , 
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where  is  an  m  x  p  matrix  and  U  -  (U^,up\  The  condition  P(X-O)  *  0 

implies  P(R2  >  0)  -1,  i.e.,  R2  satisfies  the  condition  (2.6).  As  R2  is 
independent  of  U,  by  Lemma  1  of  Anderson  and  Fang  (1982)  we  have 

A'U'DUA  -  A'UJUjA  . 

We  can  remove  A'  and  A  from  both  sides  because  Z  >  0  and  A  is  non¬ 
singular  matrix;  hence 

U'DU  3  U^UX  . 

Let  Y  be  defined  by  (4.3),  then  the  above  formula  becomes 

Y'DY/tr  Y'Y  =  Y^/tr  YY  , 

where  Y  is  partitioned  into  Y1  and  the  same  fashion  as 

and  U2 •  As  tr  Y'Y  is  independent  of  Y/tr  Y'Y,  we  can  multiply  by 
tr  Y'Y  on  both  sides  (cf.  Lemma  1  of  Anderson  and  Fang  (1982))  and  obtain 

•V 

Y'DY  3  y'Y,  'V  W  (E,m)  . 

From  Cochran's  Theorem  for  the  multivariate  normal  distribution  (cf. 

2 

Anderson  (1958)),  we  have  D-D  and  rk(D)  -  m. 

2 

Assume  D-D  and  rkD  -  m;  there  exists  an  orthogonal  matrix  T 

fl  ® ) 

such  that  I” DT  -  ~m  ~  .  Let  Z  -  TX,  then  Z  ^  LEC  (0,£,4>)  (cf. 

-  (0  0J  ~  *  nXp  ~  - 

Corollary  2  of  Lemma  2) ,  and 

x'dx  -  z'r'orz  -  z' 

. U 

Q.E.  Di 


3?- 


Hh 


*p,2(5*  V  ♦>• 
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There  Is  a  close  relationship  between  Cochran's  Theorem  in  the  uni¬ 


variate  case  and  one  In  the  multivariate  case.  When  the  population  is 
normal  this  relationship  has  been  established  (cf.  Anderson  (1958),  Kao 
(1973)).  We  will  point  out  similar  results  for  elliptically  contoured 
distributions . 

Corollary  1.  Assume  X  and  D  satisfy  the  condition  of  Theorem  9. 

Then  X'DX  %  MGp  2(Z;  <J>)  if  and  only  if  rx'DXi'v  G2(|;  <J>) 

for  every  1  a  3lp  with  £'Z£  ■  1,  where  $  €  *n •*-*■  R  =  R^n/2  n(p-l)/2’ 

Proof.  If  X'DX  ^  MG  „(Z;  <f>) ,  then  D2  -  D  and  rk(D)~m 

-  -  —  p,2  ~  i  2.  ~ 

by  Theorem  9.  From  Corollary  1  of  Lemma  2  X£  ^  LEC  .(0,1'Z £,4>)  *  EC  (0,1  ,<)>) 

—  nxi  ~  ~  ~  ~  n  -  ~n 

for  all  £  €  1RP  and  £'I£  *  1,  where  $  £  9n+-*-R*  =  R*^n/2  n(p-l)/2* 

Thus  £'X'DX£  'v-  G^j;  <J>)  by  Theorem  1  of  Anderson  and  Fang  (1982). 

If  £'X'DX£  ^G2(|;  ♦)  for  some  £  £]RP  and  £'Z£  -  1 

if  2 

with  <j>  ■*-*■  R  «  Rbn/2  n(p-l)/2'  t^en  D  “  D  and  rkH  ■  m  by  Theorem 
1  of  Anderson  and  Fang  (1982).  The  assertion  follows  by  Theorem  9.  Q.E.D. 

By  using  a  technique  similar  to  the  one  for  proving  Theorem  9,  we 
can  obtain  the  following  theorem. 

Theorem  10.  Assume  that  X  LEC _ (0,Z,$)  with  Z  >  0  and 

" "  ”  ~  nxp  <**  *•  ^ 

P(X«0)  -  0,  and  D^,...,Dk  are  symmetric  matrices;  then 

(X'DjX . X'DjX)  ^  MGpfk+i(2;  V2 . nk/2;nk+i/2;$)  where  n±  >  p, 

lc+l  2 

i  -  l,...,k,  n^+1  _>  1,  Z*  ni"  n*  **  and  °®ly  if  Di  “  "  ni* 

i  ■  1 . k  and  D.D.  -  0  for  i  t*  j. 

-i-J 
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7.  Applications. 

In  this  section  we  apply  the  theory  of  multivariate  elllptlcally 
contoured  distributions  to  the  multiple  regression  model. 

We  consider  the  following  model 


ft 


(7.1) 


«  X  B  +  E 
n*p  ~nxq~qxp  ~nxp 


(E(VecE)  *  0  ,  D(vecE') 


Minimizing  |e'E|  or  tr  E'E 
squares 


p+q  <_  n  ,  rkX  ”  q 


I  ©E,  E  ^  LEC  (0,E ,d>)  with  E  >  0. 
~n  w  ~  «.  nxp  \  ~ 


with  respect  to  B  gives  the  least 


(7.2) 


B  -  (X'X^X’Y  , 


since 


(7.3)  E'E  -  (Y-XB) ' (Y-XB)  -  (Y-XB) ' (Y-XB)  +  (B-B) 'X'X(B-B)  . 
By  assumption  we  have 

(7.4)  E  -  RUA  , 
where  A' A  “  E.  Thus 

(7.5)  Y  -  XB  +  RUA 

and 

(7.6)  B  ^  B  +  R(X'X)_1X'UA  . 
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Theorem  11.  Under  Che  assumption  in  the  model  (7.1)  we  have 

(7.7)  vec(B-B)  ^  EC  (0,1  ©  (X'X)"1,^)  , 

-  pq -  -  ~ 

where  <J>  €  *  R*  -  Rb  .  #— . 

nq  pq/2, (n-q)p/2 

Proof .  From  (7.4)  and  (7.6)  we  have  B-B  =  (X'X)~^X'E.  Then 
Theorem  11  follows  by  Corollary  2  of  Lemma  2.  Q.E.D. 

A 

In  general  the  distribution  of  B  does  not  belong  to  the  class  of 

-1  2 

the  multivariate  elliptically  contoured  distributions  unless  (X'X)  *  c  I 

where  c  is  a  constant  in  which  case 

(7.8)  B-B  ^  LEC  (0,c2Z,<J>)  , 

~  QaU  -s,  «V 

£ 

where  <p  €  S^^R  being  defined  by  Theorem  11. 


Corollary .  Denote  B  =  (8^, . . . ,8^) .  Under  the  assumption  of  Theorem 
11  we  have: 

«)  gj  •>-  EC^.a^U’X)-1.*),  whet.  ♦  6  «,"***  -  %2>npHl/2- 
(2)  cov(8i28j)  -  const.  oij(X,X)  *  . 

Proof.  The  assertion  (1)  is  a  consequence  of  Lemma  2  and  Theorem  11. 
The  assertion  (2)  follows  by  Theorem  4  of  Cambanis,  Huang  and  Simons 
(1981).  Q.E.D. 
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Now  we  consider  the  distribution  of  E'E  =  (Y-XB) ' (Y-XB) .  As 

E'E  -  Y'(I  -X(X'X)_1X')Y  -E'(I  -X(X'X)“1X') E 

~  ~  ~  ~n  ~  ~  —  ~  ~  ~n  ~  ~  ~  ~  ~ 

and  (In-X(X’X)  ^X')  is  a  projection  matrix  with  rank  n-q,  if 
P(E*0)  “  0,  then  from  Theorem  9, 

(7.9)  S  H  E'E  ^  MG  ~(I ; (n-q)/2;q/2 ;$)  . 

~  ~  P»  ~~ 

Under  the  model  (7.1)  we  want  to  test  the  following  linear  hypothesis: 

(7.10)  H:  HB  ■  C  ,  H:  txq  ,  C:  txp  and  rkH  ■  t <  p  . 

Under  the  condition  of  HB  ■  C  minimizing  |e'e|  or  tr  E'E  with 
respect  to  B  gives  the  least  squares  estimator 

(7.11)  Br  -  B  -  (X'X)-1H'(H(X'X)-1H,)_1(HB-C)  , 

A 

where  B  is  given  by  (7.2).  Since 

(Y-XB) ' (Y-XB)  -  (Y-xS  ) ' (Y-XB  )  +  (B  -B) 'X'X(BU-B)  . 

**  *"*  ^  ^  — n  *v  -H  — -  —  —  — n  ** 

Thus 

min  | (Y-XB) ' (Y-XB) |  -  j (Y-XBU) ' (Y-XB„) | 

HB-C  -  ~  ~~  ~  - H 

-  | (Y-XB) '(Y-XB)  +  (HB-C) '(H(X’X)"1H')"1(HB-C)| 

=  |S  +  T|  (say). 
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A  statistic  for  testing  the  hypothesis  H  is 


(7.12) 


X 


min| (Y-XB) ' (Y-XB) | 

B  -  ~  ~~ 


Noting  E  ”  RUA  and  B  =  B  +  (X'X)  ^X'UA,  if  the  hypothesis  is  true  we 
have 

s  =  R2A'U' [I-X(X'X)“1X']UA 

and 

T-  (HB-C) ,(H(X*X)-1H,)~1(HB-C) 

-  (B-B) ,H'(H(X'X)~1H,)~1H(B-B) 

-  R2A,U'X(X,X)“1H'(H(X,X)“1H,)-1H(X'X)-1X,UA  , 


where  A'A  “  Z  (cf.  (7.6)). 

•v  <v  «*» 

2 

Substituting  them  into  (7.12)  we  see  that  X  is  independent  of  R  . 

By  the  method  applied  in  Section  A,  the  distribution  of  X  is  the  same 

as  in  the  normal  case;  i.e.  U  .  (Wilks'  distribution). 

p,t ,n-q 


Theorem  12.  Under  the  model  (7.1)  and  P(E-O)  -  0,  the  statistic 

X  given  by  (7.12)  for  testing  the  linear  hypothesis  (7.10)  is  distributed 

according  to  U  .  _  _  . 

p,t,n-q 


References 


Anderson,  T.  W.  (1958) ,  An  Introduction  to  Multivariate  Statistical 
Analysis.  John  Wiley  and  Sons,  New  York. 

Anderson,  T.  W.,  and  Fang,  K.  T.  (1982),  Distributions  of  quadratic  forms 
and  Cochran's  Theorem  for  elliptically  contoured  distributions  and 
their  applications.  Technical  Report  No.  53,  ONR  Contract  N00014-75- 
C-0442,  Department  of  Statistics,  Stanford  University,  Stanford,  CA. 

Cambanls,  S.,  Huang,  S.,  and  Simons,  G.  (1981),  On  the  theory  of  ellip¬ 
tically  contoured  distributions,  J.  Multivariate  Analysis.  11, 

368-385. 

Devlin,  S.  J.,  Gnanadesikan,  R. ,  and  Kettenring,  J.  R.  (1976),  Some 
multivariate  applications  of  elliptical  distributions.  In  Essays 
in  Probability  and  Statistics  (S.  Ideka  et  al.,  Eds.),  365-395, 

Shink  Tsusho  Co.,  Ltd.  Tokyo. 

Johnson,  N.  J.,  and  Kotz,  S.  (1972),  Distributions  in  Statistics:  Continuous 
Multivariate  Distributions,  John  Wiley  and  Sons,  New  York. 

Kariya,  T.,  and  Eaton,  M.  L.  (1977),  Robust  tests  for  spherical  symmetry, 
Ann,  of  Statistics,  J>,  206-215. 

Kelker,  D.  (1970),  Distribution  theory  of  spherical  distributions  and  a 
location  scale  parameter  generalization,  Sankhya  A,  32,  419-430. 

Muir head ,  R.  J.  (1980),  The  effects  of  elliptical  distributions  on  some 
standard  procedures  involving  correlation  coefficients:  A  review. 
Multivariate  Statistical  Analysis  (R.  P.  Gupta  ed.),  143-159. 

Rao,  C.  R.  (1973),  Linear  Statistical  Inference  and  Its  Applications, 

John  Wiley  and  Sons,  New  York. 

Schoenberg,  I.  J.  (1938),  Metric  spaces  and  completely  monotone  functions, 
Ann.  Math..  39,  811-841. 


42 


TECHNICAL  REPORTS 


OFFICE  OF  NAVAL  RESEARCH  CONTRACT  N0001U-67-A-0112-0030  (NR-0U2-03U) 


1.  "Confidence  Limits  for  the  Expected  Value  of  an  Arbitrary  Bounded  Random 
Variable  with  a  Continuous  Distribution  Function,"  T.  W.  Anderson, 

October  1,  1969* 

2.  "Efficient  Estimation  of  Regression  Coefficients  in  Time  Series,"  T.  W. 
Anderson,  October  1,  1970. 

3.  "Determining  the  Appropriate  Sample  Size  for  Confidence  Limits  for  a 
Proportion,"  T.  W.  Anderson  and  H.  Burstein,  October  15,  1970. 

1».  "Some  General  Results  on  Time-Ordered  Classification,"  D.  V.  Hinkley, 

July  30,  1971. 

5.  "Tests  for  Randomness  of  Directions  against  Equatorial  and  Bimodal 
Alternatives,"  T.  W.  Anderson  and  M.  A.  Stephens,  August  30,  1971. 

6.  "Estimation  of  Covariance  Matrices  with  Linear  Structure  and  Moving 
Average  Processes  of  Finite  Order,"  T.  W.  Anderson,  October  29,  1971. 

7.  "The  Stationarity  of  an  Estimated  Autoregressive  Process,"  T.  W. 

Anderson,  November  15,  1971. 

8.  "On  the  Inverse  of  Some  Covariance  Matrices  of  Toeplitz  Type,"  Raul 
Pedro  Mentz,  July  12,  1972. 

9.  "An  Asymptotic  Expansion  of  the  Distribution  of  "Studentized"  Class¬ 
ification  Statistics,"  T.  W.  Anderson,  September  10,  1972. 

10.  "Asymptotic  Evaluation  of  the  Probabilities  of  Misclassification  by 
Linear  Discriminant  Functions,"  T.  W.  Anderson,  September  28,  1972. 

11.  "Population  Mixing  Models  and  Clustering  Algorithms,"  Stanley  L. 

Sclove,  February  1,  1973. 

12.  "Asymptotic  Properties  and  Confutation  of  Maximum  Likelihood  Estimates 
in  the  Mixed  Model  of  the  Analysis  of  Variance,"  John  James  Miller, 
November  21,  1973- 

13.  "Maximum  Likelihood  Estimation  in  the  Birth-and-Death  Process,"  Niels 
Keiding,  November  28,  1973. 

lU.  "Random  Orthogonal  Set  Functions  and  Stochastic  Models  for  the  Gravity 
Potential  of  the  Earth,"  Steffen  L.  Lauritzen,  December  27,  1973. 

15.  "Maximum  Likelihood  Estimation  of  Parameter'  of  an  Autoregressive 
Process  with  Moving  Average  Residuals  and  Other  Covariance'  Matrices 
with  Linear  Structure,"  T.  W.  Anderson,  December,  1973. 

16.  "Note  on  a  Case-Study  in  Box^Jenkins  Seasonal  Forecasting  of  Time  series, 
8teffen  L.  Lauritzen,  April,  197^. 


TECHNICAL  REPORTS  (continued) 


17.  "General  Exponential  Models  for  Discrete  Observations," 

Steffen  L.  Lauritzen,  May,  1974. 

13.  "On  the  Interrelationships  among  Sufficiency,  Total  Sufficiency  and 
Some  Related  Concepts,"  Steffen  L.  Lauritzen,  June,  1974. 

19.  "Statistical  Inference  for  Multiply  Truncated  Power  Series  Distributions," 
T.  Cacoullos,  September  30,  1974. 


Office  of  Naval  Research  Contract  N00014-75-C-0442  (NR-042-034) 


20.  "Estimation  by  Maximum  Likelihood  in  Autoregressive  Moving  Average  Models 
in  the  Time  and  Frequency  Domains,"  T.  W.  Anderson,  June  1975. 

21.  "Asymptotic  Properties  of  Some  Estimators  in  Moving  Average  Models," 

Raul  Pedro  Mentz,  September  8,  1975. 

22.  "On  a  Spectral  Estimate  Obtained  by  an  Autoregressive  Model  Fitting," 

Mituaki  Huzii,  February  1976. 

23.  "Estimating  Means  when  Some  Observations  are  Classified  by  Linear 
Discriminant  Function,"  Chien-Pai  Han,  April  1976. 

2k.  "Panels  and  Time  Series  Analysis:  Markov  Chains  and  Autoregressive 
Processes,"  T.  W.  Anderson,  July  1976. 

25.  "Repeated  Measurements  on  Autoregressive  Processes,"  T.  W.  Anderson, 

September  1976. 

26.  "The  Recurrence  Classification  of  Risk  and  Storage  Processes," 

J.  Michael  Harrison  and  Sidney  I.  Resnick,  September  1976. 

27.  "The  Generalized  Variance  of  a  Stationary  Autoregressive  Process," 

T.  W.  Anderson  and  Raul  P. Mentz,  October  1976. 

28.  "Estimation  of  the  Parameters  of  Finite  Location  and  Scale  Mixtures," 

Javad  Behboodian,  October  1976. 

29.  "Identification  of  Parameters  by  the  Distribution  of  a  Maximum 

Random  Variable,"  T.  W.  Anderson  and  S.G.  Ghurye,  November  1976. 

30.  "Discrimination  Between  Stationary  Guassian  Processes,  Large  Sample 

Results,"  Will  Gersch,  January  1977. 

31.  "Principal  Components  in  the  Nonnormal  Case:  The  Test  for  Sphericity," 

Christine  M.  Watemaux,  October  1977. 

32.  "Nnnnegatlve  Definiteness  of  the  Estimated  Dispersion  Matrix  in  a 

Multivariate  Model,"  F.  Pukelsheim  and  George  P.H.  Styan,  May  1978. 


TECHNICAL  REPORTS  (continued) 


33.  "Canonical  Correlations  with  Respect  to  a  Complex  Structure," 

Steen  A.  Andersson,  July  1978. 

34.  "An  Extremal  Problem  for  Positive  Definite  Matrices,"  T.U.  Anderson  and 

I.  Olkin,  July  1978. 

35.  "  Maximum  likelihood  Estimation  for  Vector  Autoregressive  Moving 
Average  Models,"  T.  U.  Anderson,  July  1978. 

36.  "Maximum  likelihood  Estimation  of  the  Covariances  of  the  Vector  Moving 

Average  Models  in  the  Time  and  Frequency  Domains,"  F.  Ahrabi,  August  1978. 

37.  "Efficient  Estimation  of  a  Model  with  an  Autoregressive  Signal  with 
White  Noise,"  ¥.  Hosoya,  March  1979. 

38.  "Maximum  Likelihood  Estimation  of  the  Parameters  of  a  Multivariate 

Normal  Distribution, "T.W.  Anderson  and  I.  Olkin,  July  1979. 

39.  "Maximum  Likelihood  Estimation  of  the  Autoregressive  Coefficients  and 

Moving  Average  Covariances  of  Vector  Autoregressive  Moving  Average  Models," 
Fereydoon  Ahrabi,  August  1979. 

40.  "Smoothness  Priors  and  the  Distributed  Lag  Estimator,"  Hirotugu  Aka ike, 

August,  1979. 

41.  "Approximating  Conditional  Moments  of  the  Multivariate  Normal  Distribution," 

Joseph  G.  Deken,  December  1979.  ' 

42.  "Methods  and  Applications  of  Time  Series  Analysis  -  Part  Is  Regression, 

Trends,  Smoothing,  and  Differencing,"  T.W.  Anderson  and  N.D.  Singpurwalla, 
July  1980. 

43.  "Cochran's  Theorem,  Rank  Additivity,  and  Tripotent  Matrices."  T.W.  Anderson 

and  George  P.H.  Styan,  August,  i960. 

44.  "On  Generalizations  of  Cochran's  Theorem  and  Projection  Matrices," 

Akimichi  Takemura,  August,  1980. 

45.  "Existence  of  Maximum  Likelihood  Estimators  in  Autoregressive  and 

Moving  Average  Models,"  T.W.  Anderson  and  Raul  P.  Mentz,  Oct.  1980. 

46.  "Generalized  Correlations  in  the  Singular  Cace,"  Ashis  Sen  Gupta, 

November  1980. 

47.  "Updating  a  Discriminant  Function  on  the  Basis  of  Unclassified  Data," 

G.J.  McLachlan  and  S.  Ganesalingam,  November  1980. 

48.  "A  New  Proof  of  Admissibility  of  Tests  in  the  Multivariate  Analysis 

of  Variance."  T.W.  Anderson  and  Akimichi  Takemura,  January,  1981. 

49  "A-Optimality  for  Regression  Designs,"  N.  N.  Chan,  January  1981. 

50.  "Tests  for  Standardized  Generalized  Variances  of  Multivariate  Lormal 

Populations  of  Possibly  Different  Dimensions,"  Ashis  Sea  Gup~a, 

February  1981. 

51.  "A  Linear  Theory  for  Noncausality,"  J.P.  Florens  and  M.  houchart, 

November  1981. 

52.  "On  the  Problems  of  Construction  and  Statistical  Inference  Associated 

With  a  Generalization  of  Canonical  Variables,"  Ashis  Sen  Gupta, 

February  1982. 


TECHNICAL  REPORTS  (continued) 


53.  "Distributions  of  Quadratic  Forms  and  Cochran's  Theorem  for  Elliptically 
Contoured  Distributions  and  Their  Applications,"  T.  W.  Anderson  and 
Kai-Tai  Fang,  May  1982. 

54.  "On  the  Theory  of  Multivariate  Elliptically  Contoured  Distributions  and 
Their  Applications,"  T.  W.  Anderson  and  Kai-Tai  Fang,  May  1982. 

55.  "On  Tests  for  Equicorrelation  Coefficient  and  the  Generalized  Variance  of 
a  Standard  Symmetric  Multivariate  Normal  Distribution,"  Ashis  Sen  Gupta, 
May  1982. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAOE  flffMa  Oat*  Eafam 0 


REPORT  DOCUMENTATION  PAGE 


t  w.  rrca r  y*?4  /  mm 


4.  TITLE  ran*  SubtltU > 


On  the  Theory  of  Multivariate  Elliptlcally  Con¬ 
toured  Distributions  and  their  Applications 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


S.  RECIPIENT'S  CATALOG  NUMBER 


s.  tyre  of  report  a  period  covered 


Technical  Report 


S.  PERFORMING  ORO.  REPORT  NUMBER 


T.  AUTHORS 


T.  U.  Anderson  and  Kai-Tai  Fang 


ONTRACT  OR  GRANT  NUMBER?*; 


N00014-75-C-0442 


S.  PERFORMING  ORGANIZATION  NAME  ANO  AOORESS 

Department  of  Statistics 
Stanford  University 
Stanford,  California 


II.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Office  of  Naval  Research 

Statistics  and  Probability  Program  Code  436 
Arlington,  Virginia  22217 


I 


12.  REPORT  DATE 

May  1982 


IS.  NUMBER  OF  PAGES 

42 


MONITORING  AGENCY  NAME  A  ADDRESS?!'  dlllmrmn I  fraa i  Controtttnt  Olll cm)  IS.  SECURITY  CLASS,  (ol  ihla  raportj 

UNCLASSIFIED 


IS.  DISTRIBUTION  STATEMENT  (at  t Mm  Rap  or!) 


Approved  for  public  release;  distribution  unlimited. 


t7.  DISTRIBUTION  ST  ATCMENT  ( el  the  mb  a  tr  met  enteaed  In  Block  SO,  It  dlttoront  hem  Mooort) 


ft.  KEY  BOROS  (Continue  on  reeeree  cldc  It  nooooomy  M  Identify  by  block 

Cochran's  theorem,  correlation  coefficients,  generalized  variance,  Hotelling's 

T2,  multiple  correlation  coefficients,  multiple  regression  model,  multivariate 
elliptlcally  contoured  distributions,  sample  covariance  matrix. 


I.  ABSTRACT  (Continue  on  roeoeoo  aide  It  mBtiwy  end  Identity  by  block  neemboa) 

In  this  paper,  the  multivariate  elliptlcally  contoured  distributions 
which  generalize  the  elliptlcally  contouisd  distribution  to  the  case  of  a 
matrix  are  defined  and  a  special  class  of  multivariate  elliptlcally  contoured 
distributions  is  studied  In  detail.  For  this  class  we  obtain  the  distribu¬ 
tions  of  the  following  statistics:  correlation  coefficients,  multiple  corre¬ 
lation  coefficients.  Hotelling's  T2,  sample  covariance  matrix,  generalized 
variance,  characr eristic  roots  of  the  covariance  matrix,  quadratic  forms,  etc. 


T  m  .IM4W  I#  »  +  *  *  CO  )4'')|*  rn  K  B'M  IHM  I  r- 


I  JANETS  1473 


EOITION  OF  I  NOV  I*  IB  OBSOLETE 
»/n  0 102*014* AMI  i 


UNCLASSIFIE 


